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Abstract 

Background: Symbiotic relationships have contributed to major evolutionary innovations, the maintenance of 
fundamental ecosystem functions, and the generation and maintenance of biodiversity. However, the exact nature of 
host/symbiont associations, which has important consequences for their dynamics, is often poorly known due to 
limited understanding of symbiont taxonomy and species diversity. Among classical symbioses, figs and their 
pollinating wasps constitute a highly diverse keystone resource in tropical forest and savannah environments. 
Historically, they were considered to exemplify extreme reciprocal partner specificity (one-to-one host-symbiont species 
relationships), but recent work has revealed several more complex cases. However, there is a striking lack of studies with 
the specific aims of assessing symbiont diversity and how this varies across the geographic range of the host. 

Results: Here, we use molecular methods to investigate cryptic diversity in the pollinating wasps of a widespread 
Australian fig species. Standard barcoding genes and methods were not conclusive, but incorporation of phylogenetic 
analyses and a recently developed nuclear barcoding gene (ITS2), gave strong support for five pollinator species. Each 
pollinator species was most common in a different geographic region, emphasising the importance of wide 
geographic sampling to uncover diversity, and the scope for divergence in coevolutionary trajectories across the host 
plant range. In addition, most regions had multiple coexisting pollinators, raising the question of how they coexist in 
apparently similar or identical resource niches. 

Conclusion: Our study offers a striking example of extreme deviation from reciprocal partner specificity over 
the full geographical range of a fig-wasp system. It also suggests that superficially identical species may be 
able to co-exist in a mutualistic setting albeit at different frequencies in relation to their fig host's range. We 
show that comprehensive sampling and molecular taxonomic techniques may be required to uncover the 
true structure of cryptic biodiversity underpinning intimate ecological interactions. 
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Background 

Symbiosis between disparate organisms has been respon- 
sible for key innovations in evolutionary history, such as 
the origin of plants [1]. It is also pivotal in several vital 
ecosystem functions, such as nitrogen fixation [2] and pol- 
lination e.g. [3], and has been implicated as a major driver 
in the generation of biodiversity [4,5]. The intimate nature 
and high partner specificity of many symbioses establishes 
conditions for strong interactions between the species in- 
volved, in terms of both population dynamics and coevolu- 
tion [6,7]. A basic but important question is "how many 
species are involved in a host/symbiont interaction?", be- 
cause species interaction strength, a major determinant of 
ecological and evolutionary dynamics, will differ between a 
system with one host and one symbiont species and one 
with multiple host or symbiont species [8]. In addition, 
when multiple symbionts compete to utilise host resources, 
this may select for more selfish or virulent behaviour. 

The obligate mutualism between fig trees and their 
pollinating wasps is a classic and much-studied example 
of symbiosis [9-11]. Fig trees are important components 
of rainforest and savannah ecosystems and provide food 
for many vertebrate and invertebrate animals [12]. Figs 
can only be pollinated by host-specific wasp species, 
which, in turn, can only reproduce by laying their eggs in 
fig flowers, upon which their developing larvae feed. This 
association was long considered to be a textbook example 
of extreme reciprocal partner specificity, with each fig 
species having a unique pollinator species [9]. However, 
this picture has been eroded by several reports of a single 
fig species hosting multiple pollinator species and a few re- 
ports of two fig species sharing pollinator species [13-21]. 

It is now clear that the paradigm of one-to-one recip- 
rocal partner specificity is no longer tenable [9,22], but 
true patterns of partner species associations remain to 
be revealed. Recent reviews have posited that perhaps a 
third or half of the >750 fig species worldwide may have 
multiple pollinator species [9,23]. However, these esti- 
mates are best regarded as informed guesses for two 
main reasons. First, most of the evidence for multiple 
pollinators is the by-product of studies conducted for 
other reasons, such as exploring wasp reproductive be- 
haviour [15], or phylogeographic histories [18]. Thus, 
these lines of evidence do not come from targeted stud- 
ies aiming to document the numbers of species involved 
in the plant/pollinator symbiosis. In particular, there is a 
striking lack of studies involving wide sampling of in- 
sects from across host plant geographic ranges [19,24]. 

Second, there is a "taxonomy gap" that hinders assess- 
ment of species associations [25-27]. There are only about 
150 described agaonid species, but the true number seems 
likely to exceed 1000 [28]. This means that investigation 
of the pollinators of many fig species begins with only a 
genus level identification and with no described wasp 



species already linked to the fig species in question. More- 
over, even when a described pollinator is known, genetic 
studies often reveal further pollinator species that are ei- 
ther morphologically cryptic within the one named entity, 
morphologically distinguishable but not previously sam- 
pled, or previously sampled but unrecognised within 
mixed species collections, e.g. [29]. 

Molecular data have already contributed substantially to 
rejection of the old 1:1 paradigm of fig/pollinator spe- 
cificity and should now play a key role, in tandem with 
morphological analysis, in establishing true patterns of 
pollinator diversity and variation in plant/pollinator inter- 
actions at local and regional scales [19]. Wide geographic 
sampling is crucial [30], because: a) some species may not 
occur in all parts of the host range, and b) if intraspecific 
genetic variation is underestimated by sampling few sites, 
it may be harder to identify the molecular 'barcoding gap' 
(or appropriate clades in phylogeny-based methods) be- 
tween species [31,32]. 

Another important issue is the choice of markers for 
molecular taxonomy. The animal barcoding approach de- 
veloped by CBOL uses a standard section of the mito- 
chondrial COI gene [33,34]. This has proved valuable in 
some studies of fig wasps, but in a recent phylogenetic 
study of 200 species from 19 genera [35], the success rate 
for PCR and direct sequencing was much higher with 
cytochrome b and these two linked mitochondrial markers 
tend to reveal very similar patterns [36] . Both markers can 
vary substantially within fig wasp (and other) species, 
which may hinder correct species delimitation, in particu- 
lar via over- splitting [37,38]. Consequently, it is desirable 
to also use a nuclear marker and to seek congruent species 
delimitation between the two genomes [39]. With nuclear 
markers, the challenge is to find one that both amplifies 
reliably across species and shows sufficient variation for 
discrimination between closely related species [40]. None 
of the currently used nuclear markers achieve the wide- 
spread utility of favoured mtDNA markers like COI and 
cytb. However, an internal transcribed spacer region 
(ITS2) of rDNA has recently been proposed as a useful 
nuclear barcoding marker for animals [41-44]. 

A further issue with molecular data is the choice of spe- 
cies delimitation method. Barcoding-type approaches 
use genetic distance data directly to identify a barcoding 
gap between the pairwise genetic distances found within 
and between species. These approaches often work well 
in practice, but threshold genetic distances vary across 
taxa, and well-known biological phenomena, such as 
introgression and selective sweeps due to Wolbachia 
bacteria, can confound the expected patterns. Mean- 
while, phylogeny-based methods aim to identify clades 
that are evolutionary significant units (ESUs) and in- 
voke the phylogenetic species concept [45]. These have 
a justifiable conceptual basis, but are yet to be as widely 
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used and are experiencing a period of relatively rapid 
methods development. 

In this study, we investigate the diversity of pollinator 
wasps associated with a single widespread fig species, 
Ficus rubiginosa. This plant is endemic to Australia and 
occurs widely in diverse habitats, including eucalypt scrub 
and rainforest, in a roughly 2500 km coastal belt that 
stretches from tropical northern Queensland to temperate 
southern New South Wales [46]. It belongs to the Ficus 
section Malvanthera, which is pollinated by wasp species 
in the genus Pleistodontes. Its only known pollinator spe- 
cies is P. imperialis, which was originally described by 
Saunders (1882). Following new wasp collections from 
several malvantheran fig species, a taxonomic revision of 
Pleistodontes was carried out [47]. This led to the descrip- 
tion of seven new Pleistodontes species overall, but no 
change to the conclusion that F. rubiginosa was associated 
with a single pollinator species {P. imperialis). 

Subsequent genetic work by Haine et al. [17] involved 
sampling P. imperialis from several regions with the aim of 
investigating the phylogeography of the species. Instead, 
cytb data revealed four deep clades, suggesting the presence 
of cryptic species. Nuclear sequence data were obtained 



from the D2 region of 28S rDNA, but this showed almost 
no variation between individuals in the samples. A second 
nuclear region, wingless, was also studied and again 
proved almost invariant. The authors concluded that, des- 
pite lack of resolution from the nuclear markers, the data 
supported the existence of four species within the 'P. 
imperialis complex Re-examination of insect specimens 
from each clade by an expert taxonomist (J-Y Rasplus, 
INRA) revealed that one putative species (clade two in 
[17]) could be distinguished by morphology (colour) alone, 
and another (clade 1) by relative head length, demonstrat- 
ing the value of integrating molecular and morphological 
information and iterative assessment of species boundaries 
e.g. [48]. However, two further putative species (clades 
3 and 4) remain morphologically cryptic. F. rubiginosa 
therefore has one of the highest diversities of pollinator 
species known for any fig species [18,19,49]. However, as 
in other studies, this has been revealed as a by-product of 
sampling wasps for other purposes. Here, we make a tar- 
geted study with the aim of using extensive sampling to: 

1. Compare the 'performance of mitochondrial (cytb 
and COI) and nuclear (ITS2) molecular markers, 



Species 1 
Species 2 
Species 3 
Species 4 
Species 5 




Species 1: South 
Queens land/NSW 



Species 1 North Queensland 



Species i- North Queensland 



Species 4 Forty Mile 
Scrub'Undara 



Li 

Figure 1 Schematic diagram of Bayesian cytb phylogeny of 41 5 P. imperialis individuals. Colours indicate five putative species though sub-clades 
indicate up to 1 1 ESUs. Posterior node probabilities are indicated for the five species (see Supplementary Information for annotated phylogeny). 
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and different species delimitation methods (distance 
and phylogeny-based) for identifying species. 

2. Uncover the full diversity of pollinator species 
associated with F. rubiginosa and the impact of 
sampling effort on their detection. 

3. Explore geographic turnover (beta diversity) of 
pollinator species across the host plant's wide 
natural latitudinal range. 

Results 

Bayesian phylogenetic analysis of 415 Pleistodontes imper- 
ialis cytb sequences identified five species when following 
the circumscription of Haine et al. [17] (Figure 1 and 
Additional file 1: Figure SI). These are the four hypothe- 
sised species found by Haine et al. [17] and one new spe- 
cies. Support values (p = 1) are high for all five of these 
main clades. However, a case can be made for up to 11 sig- 
nificant clades (Figure 1), because there are distinct sub- 
clades within the main clades. Species 1, 2 and 3 each 
exhibit two sub-clades and species 4 exhibits three. The 
different geographical distributions of these sub-clades 
(with northern and southern sub-clades in three of the 
species - see below) further suggest a hypothesis of add- 
itional taxa, above and beyond those identified by Haine 
et al. [17], as independently evolving, non-introgressing 
lineages. 

Under a five pollinator species hypothesis, visual investi- 
gation of pairwise distances (Figure 2) does not reveal a 
clear barcode gap in cytb data as intraspecific divergences 
(0-7.2%) overlap slightly with interspecific divergences 
(4.3-18.3%). However, analysis of cytb data using jMOTU 
does favour five species (Figure 3) with a barcode gap 




1 



0.00 



005 010 015 

cytb TrN+l+G pairwise distance 



0.20 



Figure 2 Modelled TrN + 1 + G pairwise distance distribution for 
415 P. imperialis cytb sequences. No barcode gap is evident. 
Intraspecific distances range between 0-7.2%; interspecific distances 
between 4.3-1 8.3%. 




MOTU discrimination in base pairs 

Figure 3 Output from jMOTU suggests five P. imperialis MOTUs 
according to percentage sequence discrimination threshold 
employed for cytb sequences. Plateau between 1 7-30 bp 
indicates barcoding gap. 



between 17-30 base pairs discrimination. GMYC analyses 
on cytb data also indicated the existence of five species. 
The GMYC model was preferred over the null model of 
uniform branching rates (GMYC logL = 304.34, null 
model logL = 297.69, p < 0.01). All of the 330 analysed 
individuals are placed in the same species by both the 
jMOTU barcoding and GMYC approaches using the cytb 
gene. 

Bayesian analysis on a subset (n = 44) of individuals 
for the 416 bp COI region suggests six major clades 
(Additional file 2: Figure S2), comprising the five major 
cytb clades but with species 1 split into northern and 
southern clades. However, support values for clades 2-5 
are all 1.00 whilst the two species 1 sub-clades have sup- 
port values of 0.53 and 0.80 indicating poor support for a 
six species hypothesis. Moreover, jMOTU analysis of COI 
indicates five species (Additional file 3: Figure S3). Visual 
investigation of pairwise distances (Additional file 4: 
Figure S4) also does not reveal a clear barcode gap in COI 
data, with intraspecific divergences of 0-5.6% overlapping 
slightly with interspecific divergences (5.0-13.1%). Sub- 
clade congruence is absolute for individuals sequenced for 
both cytb and COI. 

Bayesian phylogenetic analyses of ITS2 data also support 
a hypothesis of five species (Figure 4). However, in con- 
trast to the cytb analysis, there is a marked lack of sub- 
clade structure. Indeed, most individuals within each of 
the five species have identical sequences, while there are 
substantial differences between species, making species 
placement unequivocal using this marker. Support values 
are high for all hypothesised species nodes. Again, all 
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Figure 4 Consensus Bayesian topology from ITS2 data for 54 P. imperialis individuals. Posterior probabilities are indicated. Tip suffices 
denote major cytb clade (i.e. species) assigned to each individual by cytb analyses (e.g. CI =cytb species 1). 



species placements are congruent with those determined 
by cytb analyses. 

Figure 5 shows the frequencies of pollinators belonging 
to each of the five P. imperialis species at sites across east- 
ern Australia (see Additional file 5: Table SI for values). 
Chi-squared tests show species to be unevenly distributed 
across populations (j 2 = 558.6048, df=28, p < 2.2xl(T 16 ). 
This pattern is also evident when testing the four most 
common species individually (species 1: jf 2 = 78.1282, df = 7, 
p = 3.315xl(T 14 ; species 2: j 2 = 441, df=7, p<2.2xl(T 16 ; 
species 3: x = 203.3034, df = 7, p < 2.2xl0~ 16 ; species 4: 
/ = 121.8537, df = 7, p < 2.2xl0~ 16 ). Most sites have more 
than one (and up to four) pollinator species present, so 
there is no absolute geographical replacement of species. 
Nevertheless, three species (1, 3 and 4) each dominate in 
particular regions, whilst species 2 is co-dominant with 
species 3 in Townsville. 

Discussion 

Symbioses have been implicated in the innovation of key 
evolutionary leaps, the development of a number of funda- 
mental ecosystem functions, and as major drivers in the 
generation of biodiversity. However, the exact nature of the 
ecological relationships between symbiotic partners is often 
unknown due to the lack of targeted investigations, ta- 
xonomic impediments and insufficient sampling. Among 



classic examples of symbioses, the fig-wasp system provides 
a crucial ecological resource across tropical and savannah 
environments. Here we provide the first comprehensive tar- 
geted investigation to assess pollinator diversity across the 
geographic range of a common and widespread fig species. 
Our study supports the existence of at least five wasp spe- 
cies pollinating Ficus rubiginosa. This P. imperialis complex 
involves more pollinator species than any other fig species 
studied to date [18], but see [21,49] and our study adds to 
the growing body of evidence that the diversity of fig wasps, 
and insects in general, has been gready underestimated 
[9,23,50]. 

Our most comprehensive sampling used mitochondrial 
cytb and delineated up to 11 ESUs that might each con- 
stitute a species (Figure 1). This cytb diversity is far 
greater than previously identified and is most likely ex- 
plained by a big increase in sampling effort relative to 
that of Haine et al. [17], whose analyses highlighted four 
largely unstructured clades. Most striking among this new 
diversity are the six individuals ascribed to species 5 that 
are distinguished from other species by extended branch 
lengths in the phylogeny. These most likely represent a 
rare fifth pollinator species not sampled previously. More- 
over, we found strong sub-clade structure within some of 
the putative species and this generally corresponds to 
largely or wholly Northern versus Southern populations of 
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Figure 5 Distribution and frequencies of five Pleistodontes 

imperialis species pollinating the fig, Ficus rubiginosa, in 

eastern Australia. N denotes sample size of geographic region. 

Townsville region includes 43 morphologically identified yellow 

pollinators. The dotted line indicates the geographic range of F. 

rubiginosa modified from [46]. 
s J 

the more widespread species. Such geographic isolation is 
suggestive of restricted gene flow and could therefore indi- 
cate distinct species, but further results argue against this. 

Investigation of modelled pairwise genetic distances 
for cytb data suggests that species delimitation is not 
straightforward among these wasps due to the lack of a 
barcoding gap. However, jMOTU analyses of cytb data 
clearly supported the existence of five cryptic species 
and this algorithm identified a barcoding gap between 
17-30 base pairs. Additionally, GMYC analyses on cytb 
data give strong statistical support for the existence of 
five P. imperialis pollinators. This method adds valuable 
weight to the five species hypothesis as it statistically 
tests the crossover point between intra- and interspecific 
branch lengths derived from genetic markers. 

In order to compare cytb phylogeny with the standard 
coxl barcoding gene, a shortened COI fragment of only 
416 bp cf 658 bp; [33] was employed on a subset of 44 
individuals. The COI topology is very similar to that for 
cytb but splits species 1 into two. However, the very low 
support values (p = 0.53 and p = 0.80) offer no strong 
evidence that either of the species 1 sub-clades are more 
closely related to species 2 for this marker. Moreover, 
jMOTU analyses support five species for this marker in 
the same way as for cytb. The nuclear ITS2 phylogeny is 
similar to that for cytb in revealing the same major ESUs. 
However, sub-clade structure is absent for this marker and 
there is essentially zero variation within any of the five 



species, allowing straightforward delimitation of taxa. Inter- 
estingly, this means that ITS2 behaves far more like a bar- 
code in the wider general sense of this word than mtDNA 
markers like COI or cytb, which typically show consider- 
able within species variation even when there is a clear bar- 
coding gap. This disparity in resolution between mtDNA 
and nDNA is most likely explained by different population 
genetic properties of these markers [51], and by multigene 
concerted evolution of the repeated rDNA cluster. 

Our sampling reveals that each of the five species has a 
different geographic range (Figure 5). Whilst there is no ab- 
solute geographical replacement of species, four species ei- 
ther dominate or are codominant in particular regions. 
Species 3 is dominant in northern Queensland, while species 
4 appears dominant in central and southern Queensland 
and in the inland region around Forty Mile Scrub. The 
yellow species 2 is found only in the Townsville region, 
where it is common and co-occurs with species 3 at simi- 
lar frequencies. Finally, species 1 occupies widely disjunct 
regions in the far north of the host plant range and in 
New South Wales, where it is the only pollinator species 
recorded. Moreover, species 1 appears excluded from 
most of the intermediate regions as no individuals were 
recorded amongst 220 wasps from the Townsville and 
Forty Mile Scrub populations. Only species 5 appears to 
have no geographical stronghold, with only six individuals 
sampled overall. However, this may be an artifact of the 
low sample size for this rare species. 

From the perspective of understanding the ecological 
dynamics of symbiosis this study clearly rejects the simple 
model of reciprocally partner-specific fig and pollinator 
species. Over most of the host plant range, it is associated 
with more than one pollinator species, but the frequency 
of the different species varies considerably between re- 
gions. Given the wide latitudinal range of the host plant, 
from the wet tropics to the temperate zone, it is possible 
that pollinator diversity is related to local adaptation to 
different climates. Whatever the causes of variation in pol- 
linator identity and diversity between regions, the vari- 
ation itself has important implications for insect/ plant 
coevolution [52]. First, the association appears to be a 1:1 
match of a fig and pollinator species in the southern part 
of the range in NSW, but to typically involve two or three 
(but not always the same) co-occurring pollinator species 
at sites further North in Queensland. This could lead to 
different coevolutionary trajectories in different parts of 
the range, e.g. because multiple symbiont species may in- 
crease host/symbiont conflict as a result of antagonistic 
competition between the symbiont species [8]. 

Another set of questions arising from our improved un- 
derstanding of pollinator species diversity and distri- 
butions surrounds the local coexistence of multiple 
pollinator species utilizing the same host resources. Eco- 
logical theory suggests that this is difficult to achieve 
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without some degree of niche separation, although an in- 
triguing possibility with fig wasps is that it might occur 
through density-dependent sex ratio (and therefore popu- 
lation growth rate) variation [53]. Different sites have dif- 
ferent sets of pollinators (Figure 5), creating the possibility 
to compare different coexistence patterns in different but 
overlapping sets of species. 

There has been only limited investigation of the com- 
parative ecology of multiple sympatric pollinator species 
associated with a given fig species [13,15-18,54] and most 
such studies have not revealed clear differences. Two 
exceptions are cases where one of the pollinators is a 
"cheat" that does not actively pollinate its host plant 
[55,56], and cases where one species is a diurnal, and the 
other a nocturnal, disperser e.g. [57]. The latter case may 
be relevant in the P. imperialis complex because species 
two is yellow rather than dark brown or black, and previ- 
ous studies in Africa and Asia have shown that pale 
coloured wasps tend to be nocturnal dispersers, while 
dark ones tend to be diurnal e.g. [57]. This is one potential 
axis for niche differentiation, but four of the five P. impe- 
rials complex species are brown/black wasps so are un- 
likely to be differentiated in this way. Further research 
into host use by the P. imperialis species complex could 
focus on genetic and morphological variation within F. 
rubiginosa. 

Another key message from our study is the importance 
of a comprehensive sampling regime. Haine et al. [17] sam- 
pled 71 P. imperialis from across F. rubiginosds range and 
hypothesised the existence of four cryptic species. Our sam- 
pling covered much of the same geographic region in a 
more thorough manner, as well as some additional loca- 
tions. Despite sampling 415 wasps, we only captured six in- 
dividuals of the newly identified species 5 from three 
coastal sites up to 1000 km apart in central and southern 
Queensland. It remains possible that further 'P. imperialis 
species are yet to be sampled. The sampling effort from 
NSW is much less comprehensive than that of Queensland, 
but this so far appears justified as only a single species has 
been found from 20 specimens from several sites. 

Conclusions 

Our study offers comprehensive evidence that mutualistic 
P. imperialis wasps pollinating the fig species F. rubiginosa 
have diversified into five distinct (including some cryptic) 
species across their entire host plant range along the east 
coast of Australia. Given the absence of diagnostic mor- 
phological and ecological differences in this sister-species 
complex that may offer corroborating evidence for the as- 
sessment of species status in an integrative taxonomic 
framework, we have shown the utility of the internal tran- 
scribed spacer region (ITS2) of rDNA as a diagnostic tool 
for species delimitation among problematic taxa. Further- 
more, we have shown that four out of five of these 



identified taxa show strong patterns of regional geo- 
graphic dominance that suggests an adaptive origin to 
localised ecological conditions and offers a likely mechan- 
ism in mediating their co-existence performing the same 
ecological role in the same fig species. Moreover, where 
we identify a pattern of geographic co-dominance there is 
reason to believe that variation in pigmentation may indi- 
cate niche differentiation in diurnal/nocturnal dispersal 
patterns. In addition to this interspecific geographic struc- 
turing, our comprehensive sampling also suggests geo- 
graphic structure at the intraspecific level for three of the 
species. Thus, our results show that one-to-one specificity 
may often break down in symbioses and offer an explan- 
ation as to under what circumstances it may do so. Add- 
itionally, the increasing level of pollinator diversity at 
lower latitudes implies that coevolutionary trajectories be- 
tween host and symbiont will vary across the range. In 
summary, our understanding of the structuring of bio- 
diversity in a symbiotic context is likely to be simplistic 
and will require targeted studies using molecular taxo- 
nomic techniques and comprehensive sampling schemes 
if we are to better understand the coevolutionary dynam- 
ics underpinning intimate ecological interactions. 

Methods 

Sample collection 

Most sampling was conducted in Queensland from 2007- 
2009 along the eastern seaboard and immediate hinterland 
between Brisbane (26° 46S, 153°02E) in the south and 
Dimbulah (17° 01S, 145°19E) in the north. Sporadic sam- 
pling was also undertaken in New South Wales (NSW), 
and also from some planted trees outside their natural 
range in Victoria and South Australia between 2000 and 
2010. Near-ripe figs were collected from trees and placed 
into hatching jars with meshed lids that allowed air flow, 
prevented overheating, and prevented wasp escape. After 
48 h each fig and all its emerged wasps were placed into 
70% ethanol. Alternatively, figs were placed directly into 
alcohol and wasps were dissected out at a later date. Since 
most pollinator wasps developing in a given fig are 
siblings, we used only one wasp per syconium for DNA 
extraction to maximise the independence of samples. 
Sampling intensity was greater in Queensland than in 
New South Wales, but evidence for much lower diversity 
in New South Wales both before [17] and after this study 
(T. Sutton, pers. comm.) justifies the lower sampling effort 
in NSW. Additional file 6: Table S2 gives location, markers 
sequenced and GenBank accession numbers for all 415 
wasps used in this study. 

Molecular methods 

A Chelex method was used for DNA extraction [58] and 
two mitochondrial (cytb and COI) and one nuclear (ITS2) 
marker were amplified. NB the COI marker is a region of 
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coxl shorter than the standard barcoding region, where 
alternate primers have been designed to avoid nuclear 
pseudo-genes within certain taxa: [59] . Our widest genetic 
sampling of individuals focussed on cytochrome b (cytb) 
using mtDNA primers CB1 and CB2 [60]. A 396 bp frag- 
ment of cytb was amplified for 415 wasps in order to as- 
sign individuals to the putative species identified by Haine 
et al. [17]. Cytb was chosen as our primary marker in 
order to integrate our study with the earlier one by Haine 
et al. [17]. We also analysed a subset of individuals using 
the COI region to facilitate comparison between mtDNA 
genes and connectivity with other datasets. 

Subsets of 3-15 individuals per cytb clade were se- 
quenced for ITS2 (n = 54) [61] and a 416 bp fragment of 
COI (n = 44) (primers CI-J-1751 and CI-N-2191; see 
[62]) in order to clarify and confirm species delimitation. 
Cytb was amplified using a Techne Touchgene gradient 
machine with 3 min at 94°C, 30 cycles of 15 s at 95°C, 
20 s at 45°C, 30 s at 72°C, and a final elongation step of 
10 min at 72°C Amplification of COI consisted of 5 min 
at 94°C, 30 cycles of 30 s at 94°C, 45 s at 50°C, 60 s at 
72°C then 10 min at 72°C. Amplification of ITS2 used 
5 min at 94°C, 35 cycles of 30 s at 94°C, 40 s at 55°C, 40 s 
at 72°C then 10 min at 72°C. Subsequent purification and 
sequencing reactions were conducted by Macrogen Inc. 
Purification was performed using ethanol precipitation and 
sequencing by BigDyeTM terminator cycling conditions 
and a 3730x1 DNA analyser. All sample sequences have 
been deposited in GenBank: cytb [GenBank: KM249475- 
KM249848], COI [GenBank: KM249375 - KM249419], 
ITS2 [GenBank: KM249420 - KM249474]; see Additional 
file 6: Table S2. 

Sequence alignment and phylogenetic analysis 

Sequence quality was checked using Finch TV Version 
1.4.0 and sequences were edited and aligned using BioEdit 
[63]. In order to place collected samples into the clade 
(hereafter species) categorisations of Haine et al. [17], and 
as a template with which to align newly generated se- 
quences, 23 P. imperialis cytb sequences with geograph- 
ical location data were downloaded from GenBank/EMBI 
(from accession numbers AJ298439 and AY567594 - 
AY567638). ITS2 sequences were checked for microsatel- 
lites [64] and a small region of 10-16 bp 'TC repeats was 
identified. Analyses were run with and without this region, 
which caused no differences in phylogenetic inference, 
other than a 2% and 3% reduction in support values at 
two internal nodes. P. imperialis ITS2 sequence lengths 
were between 312-317 bp due to the presence of indels. 
Indels ranged between 0-4 bp within individual species 
with no indels found in species 2. Within species 1, a 2 bp 
indel was diagnostic of specimens from Queensland and 
New South Wales clades. No obvious pattern was evident 
among the distribution of indels for other species. 



Additionally, ITS2 sequence length was 474 bp in the out- 
group Pleistodontes sp from F. glandifera and created 
indels throughout the alignment. No evidence of pseudo- 
genes or heteroplasmy was noted for mitochondrial re- 
gions. Bayesian methods were employed to construct 
phylogenies using MrBayes [65]. The model of nucleotide 
substitution for each gene was chosen using jModelTest2 
[66] . For mitochondrial markers the TrN + I + G and 
TPMluf + G models were chosen for cytb and COI re- 
spectively whilst the TPM3uf + G model was chosen for 
ITS2. 



Species delimitation methods 

We used two main strategies, distance (barcoding-type) 
and phylogeny-based, to delimit species. We applied two 
distance methods to cytb and COI data: (i) we calculated 
pair-wise genetic distances in PAUP [67] and plotted them 
in order to determine visually whether a barcode gap was 
present; (ii) we used jMOTU software [68] to identify the 
inflection point in the frequency distribution of genetic 
distances, which is purported to show the barcoding gap. 
For jMOTU analyses, only sequences without ambiguous 
base calls were utilized. 

As a complementary phylogeny-based analysis on cytb, 
we identified evolutionary significant units (ESUs) by 
applying the statistical GMYC method of Pons et al. [69] 
to the cytb phylogeny, using the R-Splits package [70] in R 
[71]. The GMYC method requires an ultrametric phylo- 
genetic tree and this was created from cytb data with 
BEAST software [72]. Models of evolution implemented 
in BEAST were chosen using Path Sampling and Stepping 
Stone techniques [73]. These allow the calculation of 
Bayes factors, which were evaluated using the criteria of 
Raftery [74], and the HKY + G model was chosen. Two 
independent runs of 120 million generations were per- 
formed, under constant-clock conditions with a constant 
coalescent model of species evolution. Genealogies and 
model parameters were sampled once every 6000 itera- 
tions. Preliminary runs of the GMYC using only P. imper- 
ialis sequences yielded non-significant results due to a lack 
of statistical power, probably due to the limited number of 
species despite large numbers of individuals. Consequently, 
we added sequences from four other Pleistodontes species 
and this led to significant results. 

Finally, we compared the geographic distributions of the 
pollinator species across eight regional populations. We 
used Chi-square tests in R [71] to test null hypotheses that 
each species occurs at the same frequency in each region. 

Availability of supporting data 

Sequence data has been deposited on Genbank whilst 
population frequencies are in Supporting Information. 
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Additional files 



Additional file 1: Figure SI. Annotated phylogeny of cytb mtDNA 
(see Figure 1 legend for details). 

Additional file 2: Figure S2. Consensus Bayesian topology from COI 
data for 44 P. imperialis individuals. Posterior probabilities are indicated. 
Tip suffices denote major cytb clade (i.e. species) assigned to each 
individual by cytb analyses (e.g. CI =cytb species 1). 

Additional file 3: Figure S3. Figure 3. Output from jMOTU suggests five 
P. imperialis MOTUs according to percentage sequence discrimination 
threshold employed for COI sequences. Largest plateau between 18-22 bp 
indicates barcoding gap. 

Additional file 4: Figure S4. Modelled TPMluf + G pairwise distance 
distribution for 44 P. imperialis COI sequences. No barcode gap is evident. 
Intraspecific distances range between 0-5.6%; interspecific distances 
between 5.0-13.1%. 

Additional file 5: Table SI. Population frequency counts for five 
species of Pleistodontes imperialis at eight locations in eastern Australia. 

Additional file 6: Table S2. Locations, genetic markers sequenced, and 
GenBank accession numbers for 415 Pleistondotes imperialis wasps 
sequenced in the study. 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

CTD: conducted fieldwork collections, molecular genetic work, sequence 
alignments, phylogenetic inference, statistical analyses and drafted the 
manuscript. SAB: conducted fieldwork collections, molecular genetic work 
and sequence alignments. JMC: conducted fieldwork collections, drafted the 
manuscript and conceived the study. All authors read and approved the final 
manuscript. 

Acknowledgements 

We would like to thank Brent Emerson and Simon van Noort for helpful 
comments regarding this manuscript, Simon Segar and Derek Dunn for 
helping with data collection, and Jean-Yves Rasplus for taxonomic advice. 
This work was funded by a NERC studentship (#NE/G52391 2/1 - CTD). There 
was no involvement of the funding sources in the study design; in the col- 
lection, analysis and interpretation of data; in the writing of the report; and 
in the decision to submit the article for publication. 

Author details 

'School of Biological Sciences, University of Reading, Reading RG6 6AS, UK. 
department of Biology, Syracuse University, 107 College Place, Syracuse, NY 
13244, USA. 3 Royal Horticultural Society, Wisley Garden, Woking, Surrey GU23 
6QB, UK. 4 Hawkesbury Institute for the Environment, University of Western 
Sydney, Locked Bag 1797, Penrith South DC, NSW 1797, Australia. 

Received: 22 May 2014 Accepted: 13 August 2014 



References 

1. Margulis L: Origin of Eukaryotic Cells. Newhaven, Connecticut: Yale University 
Press; 1970. 

2. Kiers ET, Rousseau RA, West SA, Denison RF: Host sanctions and the 
legume-rhizobium mutualism. Nature 2003, 425:78-81. 

3. Darwin C: On the Origin of Species by Means of Natural Selection, or the 
Preservation of Favoured Races in the Struggle for Life London, UK: John 
Murray; 1859. 

4. Leigh EG: The evolution of mutualism. J Evol Biol 2010, 23(12)2507-2528. 

5. Althoff DM, Segraves KA, Johnson MTJ: Testing for coevolutionary 
diversification: linking pattern with process. Trends Ecol Evol 2014, 
29(2):82-89. 

6. Thrall PH, Hochberg ME, Burdon JJ, Bever JD: Coevolution of symbiotic 
mutualists and parasites in a community context. Trends Ecol Evol 2007, 
22(3):120-126. 



7. Rooney N, McCann KS: Integrating food web diversity, structure and 
stability. Trends Ecol Evol 2012, 27(l):40-46. 

8. Herre EA, Knowlton N, Mueller UG, Rehner SA: The evolution of 
mutualisms: exploring the paths between conflict and cooperation. 
Trends Ecol Evol 1999, 14(2)49-53. 

9. Cook JM, Rasplus JY: Mutualists with attitude: coevolving fig wasps and 
figs. Trends Ecol Evol 2003, 18(7)241-248. 

10. Herre EA, Jander KC, Machado CA: Evolutionary ecology of figs and their 
associates: recent progress and outstanding puzzles. Annu Rev Ecol Evol 
Syst 2008, 39:439-458. 

1 1 . Segar ST, Cook JM: The dominant exploiters of the fig/pollinator 
mutualism vary across continents, but their costs fall consistently on the 
male reproductive function of figs. Ecol Entomol 2012, 37:342-349. 

12. Kissling WD, Rahbek C, Bohning-Gaese K: Food plant diversity as 
broad-scale determinant of avian frugivore richness. Proc R Soc B-Biol Sci 
2007, 274(161 1):799-808. 

13. Lopez-Vaamonde C, Rasplus JY, Weiblen GD, Cook JM: Molecular 
phylogenies of fig wasps: partial cocladogenesis of pollinators and 
parasites. Mol Phylogenet Evol 2001 ,21(1 ):55-71 . 

14. Weiblen GD: How to be a fig wasp. Annu Rev Entomol 2002, 47:299-330. 

15. Molbo D, Machado CA, Sevenster JG, Keller L, Herre EA: Cryptic species of 
fig-pollinating wasps: implications for the evolution of the fig-wasp 
mutualism, sex allocation, and precision of adaptation. Proc Natl Acad Sci 
USA 2003, l00(10):5867-5872. 

1 6. Molbo D, Machado CA, Herre EA, Keller L: Inbreeding and population structure 
in two pairs of cryptic fig wasp species. Mol Ecol 2004, 1 3(6):1 61 3-1 623. 

1 7. Haine ER, Martin J, Cook JM: Deep mtDNA divergences indicate cryptic 
species in a fig-pollinating wasp. BMC Evol Biol 2006, 6:83. 

18. Chen Y, Compton SG, Liu M, Chen X-Y: Fig trees at the northern limit of 
their range: the distributions of cryptic pollinators indicate multiple 
glacial refugia. Mol Ecol 2012, 21 (7):1 687-1 701 . 

19. Cornille A, Underhill JG, Cruaud A, Hossaert-McKey M, Johnson SD, Tolley 
KA, Kjellberg F, van Noort S, Proffit M: Floral volatiles, pollinator sharing 
and diversification in the fig-wasp mutualism: insights from Ficus 
natalensis, and its two wasp pollinators (South Africa). Proc Roy Soc B: Biol 
So 2012, 279:1731-1739. 

20. McLeish MJ, van Noort S: Codivergence and multiple host species use by 
fig wasp populations of the Ficus pollination mutualism. BMC Evol Biol, 
12:1-2012. 

21 . Compton SG, Grehan K, van Noort S: A fig crop pollinated by three or 
more species of agaonid fig wasps. African Entomol 2009, 17(2)215-222. 

22. Machado CA, Robbins N, Gilbert MTP, Herre EA: Critical review of host 
specificity and its coevolutionary implications in the fig/fig-wasp 
mutualism. Proc Natl Acad Sci USA 2005, 1 02:6558-6565. 

23. Cook JM, Segar ST: Speciation in fig wasps. Ecol Entomol 2010, 35:54-66. 

24. Kobmoo N, Hossaert-Mckey M, Rasplus JY, Kjellberg F: Ficus racemosa is 
pollinated by a single population of a single agaonid wasp species in 
continental South-East Asia. Mol Ecol 2010, 19:2700-2712. 

25. Kaartinen R, Stone GN, Hearn J, Lohse K, Roslin T: Revealing secret liaisons: 
DNA barcoding changes our understanding of food webs. Ecol Entomol 
2010, 35(5):623-638. 

26. Smith MA, Eveleigh ES, McCann KS, Merilo MT, McCarthy PC, Van Rooyen Kl: 
Barcoding a quantified food web: crypsis, concepts, ecology and 
hypotheses. PLoS One 201 1, 6(7): el 4424. 

27. Wheeler QD, Raven PH, Wilson EO: Taxonomy: impediment or expedient? 
Sci (New York, NY) 2004, 303(5656)285. 

28. Cook JM, West SA: Figs and fig wasps. Curr Biol 2005, 1 5(24):R978-R980. 

29. Lopez-Vaamonde C, Dixon DJ, Cook JM, Rasplus JY: Revision of the 
Australian species of pleistodontes (Hymenoptera : Agaonidae) 
fig-pollinating wasps and their host-plant associations. Zool J Linn Soc 
2002, 136(4):637-683. 

30. Vogler AP, Monaghan MT: Recent advances in DNA taxonomy. J Zool Syst 
Evol Res 2007, 45(1):1-10. 

31 . Moritz C, Cicero C: DNA barcoding: promise and pitfalls. Plos Biol 2004, 
2(10):1 529-1531. 

32. Will KW, Rubinoff D: Myth of the molecule: DNA barcodes for species 
cannot replace morphology for identification and classification. 

Cladistics- IntJ Willi Hennig Soc 2004, 20(1)47-55. 

33. Hebert PDN, Ratnasingham S, deWaard JR: "Barcoding animal life: 
cytochrome c oxidase subunit 1 divergences among closely related 
species,". Proc R Soc London Ser B-Biological Sci 2003, 270:S96-S99. 



Darwell ef al. BMC Evolutionary Biology 2014, 14:189 
http://www.biomedcentral.com/1471-2148/14/189 



Page 10 of 10 



34. Hebert PDN, Cywinska A, Ball SL, DeWaard JR: Biological identifications 
through DNA barcodes. Proc R Soc London Ser B-Biological Sci 2003, 
270(1512)313-321. 

35. Cruaud A, Jabbour-Zahab R, Genson G, Couloux A, Peng Y-Q, Da Rong Y, 
Ubaidillah R, Santinelo Pereira RA, Kjellberg F, van Noort S, Kerdelhue C, 
Rasplus J-Y: Out of Australia and back again: the world-wide historical 
biogeography of non-pollinating fig wasps (Hymenoptera: 
Sycophaginae). J Biogeogr 201 1, 38(2)209-225. 

36. Lopez-Vaamonde C, Wikstrom N, Kjer KM, Weiblen GD, Rasplus JY, Machado 
CA, Cook JM: Molecular dating and biogeography of fig-pollinating 
wasps. Mol Phylogenet Evol 2009, 52(3):71 5-726. 

37. Avise J: Phylogeography: The History And Formation Of Species. Cambridge 
Massachusetts, London England: Harvard University Press; 2000. 

38. Wiemers M, Fiedler K: Does the DNA barcoding gap exist? - a case study 
in blue butterflies (Lepidoptera: Lycaenidae). Front Zool 2007, 4:8. 

39. Hurst GDD, Jiggins FM: Problems with mitochondrial DNA as a marker in 
population, phylogeographic and phylogenetic studies: the effects of 
inherited symbionts. Proc R Soc B-Biol Sci 2005, 272(1572):! 525-1 534. 

40. Lohse K, Sharanowski B, Blaxter M, Nicholls JA, Stone GN: Developing EPIC 
markers for chalcidoid Hymenoptera from EST and genomic data. 

Mol Ecol Resour 201 1 , 1 1 (3):521-529. 

41. Coleman AW: Is there a molecular key to the level of 'biological species' 
in eukaryotes? A DNA guide. Mol Phylogenet Fvol 2009, 50(1 ):1 97-203. 

42. Li YW, Zhou X, Feng G, Hu HY, Niu LM, Hebert PDN, Huang DW: COI and 
ITS2 sequences delimit species, reveal cryptic taxa and host specificity of 
fig-associated Sycophila (Hymenoptera, Eurytomidae). Mol Ecol Resour 
2010, 10(l):31-40. 

43. Muller T, Philippi N, Dandekar T, Schultz J, Wolf M: Distinguishing species. 
Rna-a Publ Rna Soc 2007, 13(9):1469-1472. 

44. Zhou MJ, Xiao JH, Bian SN, Li YW, Niu LM, Hu HY, Wu WS, Murphy RW, 
Huang DW: Molecular approaches identify known species, reveal cryptic 
species and verify host specificity of Chinese Philotrypesis 
(Hymenoptera: Pteromalidae). Mol Ecol Resour 2012, 12(4)598-606. 

45. Eldredge N, Cracraft J: Phylogenetic Patterns And The Evolutionary 
Process Method And Theory In Comparative Biology. In Phylogenetic 
Patterns and the Evolutionary Process: Method and Theory in Comparative 
Biology Viii+349p. Edited by Eldredge N, Cracraft J. New York, NY, USA. II I us. 
Maps: Columbia University Press; 1980. 

46. Dixon DJ, Jackes BR, Bielig LM: Figuring out the figs: the ficus obliqua- 
ficus rubiginosa complex (Moraceae : Urostigma sect. Malvanthera). 
AustSyst Bot 2001, 14(1 ):1 33-1 54. 

47. Lopez-Vaamonde C, Dixon DJ, Cook JM, Rasplus JY: Revision of the 
Australian species of pleistodontes (Hymenoptera : Agaonidae) 
fig-pollinating wasps and their host-plant associations. Zool J Linn Soc 
2002, 136(4):63 7-683. 

48. De Queiroz K: Species concepts and species delimitation. Syst Biol 2007, 
56(6):879-886. 

49. Moe AM, Weiblen GD: Molecular divergence in allopatric ceratosolen 
(Agaonidae) pollinators of geographically widespread ficus (Moraceae) 
species. Ann Entomol Soc Am 2010, 103(6):1 025-1 037. 

50. Bruno JF, Cardinale BJ: Cascading effects of predator richness. Front Ecol 
Environ 2008, 6(10):539-546. 

51. Lowe Harris A, S Ashton P: Ecological Genetics. Carlton, Victoria, Australia: 
Blackwell Science limited; 2004. 

52. Thompson JN: The Geographic Mosaic of Coevolution. Chicago: The 
University of Chicago Press; 2005. 

53. Zhang DY, Lin K, Hanski I: Coexistence of cryptic species. Ecol Lett 2004, 
7(3):165-169. 

54. Weiblen GD, Bush GL: Speciation in fig pollinators and parasites. Mol Ecol 
2002, 11(8):1 573-1 578. 

55. Compton SG, Holton KC, Rashbrook VK, Vannoort S, Vincent SL, Ware AB: 
Studies of Ceratosolen galili, a non-pollinating agaonid fig wasp. 
Biotropica 1991, 23(2):1 88-1 94. 

56. Peng Y-Q, Duan Z-B, Yang D-R, Rasplus J-Y: Co-occurrence of two 
Eupristina species on Ficus altissima in Xishuangbanna, SW China. 
Symbiosis 2008, 45(1-3)3-14. 

57. Warren M, Robertson MP, Greeff JM: A comparative approach to 
understanding factors limiting abundance patterns and distributions in a 
fig tree-fig wasp mutualism. Ecography (Cop) 2010, 33(1 }: 1 48—1 58. 

58. West SA, Cook JM, Werren JH, Godfray HCJ: Wolbachia in two insect 
host-parasitoid communities. Mol Ecol 1998, 7(1 1 ):1 457-1 465. 



59. Xiao J-H, Wang N-X, Li Y-W, Murphy RW, Wan D-G, Niu L-M, Hu H-Y, Fu Y-G, 
Huang D-W: Molecular approaches to identify cryptic species and 
polymorphic species within a complex community of fig wasps. 
PLoSOne 2010, 5(1 1): el 5067. 

60. Jermiin LS, Crozier RH: The cytochrome-B region in the mitochondrial 
DNA of the ant Tetraponera rufoniger - sequence divergence in 
hymenoptera may be associated with nucleotide content. J Mol Evol 
1994, 38(3)282-294. 

61 . White T, Bruns T, Lee S, Taylor J: Amplification And Direct Seguencing Of 
Fungal Ribosomal Genes For Phylogenies. New York, N.Y.: Academic; 1990. 

62. Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P: Evolution, 
weighting, and phylogenetic utility of mitochondrial gene-sequences 
and a compilation of conserved polymerase chain-reaction primers. 
Ann Entomol Soc Am 1994, 87(6):651 -701 . 

63. Hall TA: BioEdit: a user friendly biological sequence alignment editor and 
analysis program for Windows95/98NT. Nucleic Acids Symp 1999, 41:95-98. 

64. Harris DJ, Crandall KA: Intragenomic variation within ITS1 and ITS2 of 
freshwater crayfishes (Decapoda: Cambaridae): implications for 
phylogenetic and microsatellite studies. Mol Biol Evol 2000, 17:284-291. 

65. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference 
under mixed models. Bioinformatics 2003, 1 9(1 2):1 572-1 574. 

66. Darriba D, Taboada GL, Doallo R, Posada D: jModelTest 2: more models, 
new heuristics and parallel computing. Nat Methods 2012, 9(8):772. 

67. Swofford DL: PAUP* Phylogenetic Analysis Using Parsimony ("and Other 
Methods). Sinauer Associates; Sunderland, Massachusetts, USA 2002. 

68. Jones M, Ghoorah A, Blaxter M: jMOTU and taxonerator: turning DNA 
barcode sequences into annotated operational taxonomic units. 
PLoSOne 201 1, 6(4):10. 

69. Pons J, Barraclough TG, Gomez-Zurita J, Cardoso A, Duran DP, Hazell S, 
Kamoun S, Sumlin WD, Vogler AP: Sequence-based species delimitation 
for the DNA taxonomy of undescribed insects. Syst Biol 2006, 
55(4)595-609. 

70. Ezard T, Fujisawa T, Barraclough T: Species' Limits by Threshold Statistics. 2009. 

71 . R. R Development Core Team: R: A language and environment for 
statistical computing. R Found Stat Comput 201 1, 1:409. 

72. Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by 
sampling trees. BMC Evol Biol 2007, 7:214. 

73. Baele G, Lemey P, Rambaut A, Suchard MA, Alekseyenko AV: Improving the 
accuracy of demographic and molecular clock model comparison while 
accommodating phylogenetic uncertainty. Mol Biol Evol 2012, In press. 

74. Raftery AE: Hypothesis testing and model selection. In Markov chain 
Monte Carlo in practice. Edited by Gilks WR, Richardson S, Spiegelhalter DJ. 
London: Chapman & Hall; 1996. 



doi:1 0.1 1 86/sl 2862-01 4-01 89-9 

Cite this article as: Darwell et al:. Molecular species delimitation of a 
symbiotic fig-pollinating wasp species complex reveals extreme 
deviation from reciprocal partner specificity. BMC Evolutionary Biology 
2014 14:189. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



